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" , Abstract 

. In this article, we calculate the contribution from the nonfactorizable soft hadronic 

' matrix element to the decay — > Xcitt" with the light-cone QCD sum rules. The 

. numerical results show that its contribution is rather large and should not be neglected. 

' The total amplitudes lead to a branching fraction which is in agreement with the 



(N 
> 



O 



experimental data marginally. 
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p_^• 1 Introduction 

D ■ 

' Recently, the Belle Collaboration measured the branching fraction for the Cabibbo- and 

color-suppressed decay Xci7i"° based on a data sample of 657 x 10^ BB events 

collected at the T(4S') resonance with the Belle detector at the KEKB asymmetric-energy 
1/^ ' e^e~ collider [1]. The signal is 40ib9 events with a significance of 4.7(J including systematic 

On ■ uncertainties, and the branching fraction is about (1.12 it 0.25 it 0.12) x 10^^. 

The decay takes place through the process b dec (or more precise b — > dec, they 
^ relate with each other by charge conjunction, in this article, we calculate the amplitudes 

■ for the process b — > dec, then take charge conjunction to obtain the branching fraction.) at 

OO , the quark-level [2j. If the leading order tree diagram dominates, the time-dependent CP- 

violating asymmetries are predicted to be the same as the ones in the b — > sec decays [3]. 
The time-dependent CP-violation parameters for the similar decay B^ — > J/ipir^ have been 
^ . measured by the Belle [H [5] and Babar [U El [8] Collaborations. The deviation of the CP- 

H ' violating asymmetries from those expectations may indicate non-negligible contributions 

from the penguin amplitudes or new physics. 

The quantitative understanding of the i?-decays depends on our knowledge about the 
nonperturbative hadronic matrix elements of the operators entering the effective weak 
Hamiltonian [2]. In recent years, great progresses have been made in this aspect, such as 
the generalized factorization approach [OllTO], the QCD factorization approach [111 I12j. 
the perturbative QCD [13t il4j . the soft-collinear effective theory [15j, etc. Factorization 
of the hadronic matrix elements has been proved to hold in the leading order in many 
processes. 

The effects of the soft gluons which violate factorization are supposed of order 0( ^^^ ) 
and neglected in the QCD-improved factorization studies [11^ il2j. however, no theoretical 
work has ever proved that they are small quantities. For the color-suppressed B to char- 
monia decays, there may be significant impacts of the nonfactorizable soft contributions. 

In Ref.|16j. Khodjamirian introduce a technique based on the light-cone QCD sum 
rules to estimate the nonfactorizable soft contributions, where the soft gluon effects are 
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represented by the quark-antiquark-gluon distribution amplitudes of the light mesons, the 
hadronic matrix element appears as a part of the hadronic dispersion relation for the 
correlation function. Thereafter, the light-cone QCD sum rules are applied to study the 
nonfactorizable hadronic matrix elements in the -B-decays due to the soft gluons exchanges 
[I71[ig[ia|20l[2ll|22l|23l[Sll25l[2S[27]. 

It is interesting to study the nonfactorizable soft contributions in the decay Xcit^^ 
with the light-cone QCD sum rules. 

The article is organized as follows: the factorizable contributions from the effective 
weak Hamiltonian are derived in Sec. 2; the soft hadronic matrix element (xci7r^|C'(0)|i?) is 
calculated with the light-cone sum rules approach in Sec. 3; numerical results are presented 
in Sec. 4; the section 5 is reserved for conclusion. 



2 Effective weak Hamiltonian and factorizable contributions 



The effective weak Hamiltonian for the b — > dec decay modes can be written as (for detailed 
discussion of the effective weak Hamiltonian, one can consult Ref. [2]) 

10 ^ 



^ I 1=3 



(1) 



where Vij^s are the CKM matrix elements, Cj's are the Wilson coefficients calculated at 
the renormalization scale /i ~ 0{mj,) and the relevant operators Oi are given by 



Oi 

3(5) 



{dabf3)v-Aici3Ca)v-A , 
{daba)v-Aici3Ci3)v-A , 

idaba)v-A'^{qpqp)v-A{V+A) j 



4(6) 



{dabf3)v^A "^{qpq 



'a)V~A(V+A) , 



(2) 



where we have neglected the Wilson coefficients C7, Cg, Cg, Cio due to their small values. 
We can reorganize the color-mismatched quark fields into color singlet states by Fierz 
transformation, and express the effective weak Hamiltonian in the following form. 



where 



V2 



C2{p.) + ^4r^]02 + 2Ci{ii)02 



O2 = (c7M(l-75)Yc)(d7'^(l-75)Y&), 
0, = (C7^(l-75)^c)(d7'^(l-75)y?'), 
O, = (c7M(l + 75)yc)(d7'^(l-75)y&), 



(3) 



(4) 
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and A°'s are SU{3) Gell-Mann matrices. 

The factorizable matrix elements {xciip)'^^ {—q)\Hw\B (p — q)) can be parameterized as 
e* ■ qAf, 



e* ■ qAf 



{Xci{p)\ci^,l^c\Q){7:\-q)\d^%\B{p - q)) . 



(5) 



The Xci meson decay constant is defined by (Xci(p)|c(0)7^75c(0)|0) = fxci^xa^li- '^^^ 
B — ir^ form-factor can be parameterized as 



2{7T^{q)\d^^b\BiP)) = f(p')(p + q)^-[f(p')-f,(p')]. 
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(6) 



the form-factors f{p'^) and fo{p'^) can be estimated with the hght-cone sum rules approach, 
here we take the value /("z^^J = 0.62 [28] . 

The concise expression for the factorizable matrix elements can be written as 



Ai 



-Gpfxci^xcifi^xci) 1 ^cbVc 



C2(/X) + 



VtbVt 



td 



C4(m) - Ceifl) 



(7) 



3 Light-cone QCD sum rules for the nonfactorizable hadronic 
matrix element {xci'^'^\O{0)\B) 

In the following, we apply the approach developed in Ref.[T6] for the B — > vrvr channel to 
estimate the contribution from the soft-gluon exchanges in the decay B^ — > Xcit^^- We 
write down the correlation function Ilp{p,q,k) firstly. 



Up{p,q,k) 



/ d4ye*(P-'=)^(0|T{J,(y)O(0)J5(x)}K0(<7)) , (8) 



where O = O2, the currents Jp = c^p^^c and J5 = mi^bi^^d interpolate the mesons Xci 
and B, respectively. 

The correlation function Ilp{p, q, k) can be calculated by the operator product expan- 
sion approach near the light-cone ~ ~ (x — y)^ ~ in perturbative QCD theory. It 
is function of three independent momenta chosen to he q, p — k and k. We introduce the 
unphysical momentum k in order to avoid that the B meson has the same four-momentum 
before {p — q) and after the decay (P), and thus avoid a continuum of light contributions 
in the dispersion relation in the 5-channel. The independent kinematical invariants can 
be taken as {p - q)'^, {p - fe)^, , k'^ , P'^ = {p - k - qf and p^ . We set A;^ = and take 
gi2 = _ neglecting the small corrections of order 0(m^/m^). The momentum p^ is 
kept undefined for the moment in order to make the derivation of the sum rules without 
restriction. Its value is going to be set later in this section, and chosen to be p^ = "m-^^-^ . 
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The values of {p — k)"^, {p — g)^ and should be spacelike and large in order to stay far 
away from the hadronic thresholds in the B and Xci channels. All together, we have 



k'^ = 0, / 



m 



Xc 



^, |(p-A;)|2 > AQC7D,|(p-g)|^ » Aqcd,|^'|^ > Aqcd • 



The correlation function Ilp{p, q, k) can be decomposed as 

np(p,g,A:) = {p-k)pIl + qpIli + kpU.2 + epo,p„qy^k''Ii^. 



(9) 



due to Lorentz covariance. 

According to the basic assumption of quark-hadron duality in the QCD sum rules 
|29l I30j . we insert a complete set of intermediate states with the same quantum numbers 
as the current operators J^{y) and J^i^x) into the correlation function Tlp{p, q, k) to obtain 
the hadronic representation. After isolating the pole terms of the ground state mesons Xci 
and B, we get the following result, 



{0\JpiO)\xciip - k)) 



m 



Xcl 



{p - ky^ 



{Xciip-k)\OiO)\Bip-q))\7r^iq)) 



,...{B{p-q)\MOm 



m 



B 



+ ■ 



fBm% 



{p — kY m?n — {p — q) 



m 



:{Xciip-k)7r\-q)\OmBip-q))+.^ 



Xcl 



fBm% 



m 



Xcl 



{p — kY m?n — {p 



-gap + 



{p- k)a{p- k)p 



{p - ky 



(10) 



where we have used the definition {xci{p — k)Tr^ {—q)\0 {0)\B {p — q)) = e* ■ qAn, and do 
not show the contributions from the higher resonances and continuum states above the 
corresponding thresholds explicitly, they can be written in terms of dispersion integrals 
and the spectral density can be approximated by the quark-hadron duality ansatz. 

Now we carry out the operator product expansion near the light-cone to obtain the 
representation at the level of quark-gluon degrees of freedom. Let us write down the 
propagator of a massive quark in the external gluon field in the Fock-Schwinger gauge 
firstly [31J, 



Sij{x,y) 



(27r)4 
+ m 



^ + m 



m 



2 (A;^ — m?Y k"^ — m? 



dv gs Gf^{vx + {l-v)y) 
v{x - y)p-iu 



fill 



where Ga'^ is the gluonic field strength, gg denotes the strong coupling constant. 

Substituting the above b and c quark propagators into the correlation function Hp, 
and completing the corresponding integrals, we can obtain the hadronic spectral density 
at the level of quark-gluon degrees of freedom. In calculation, the following three-particle 
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vr'^ light-cone distribution amplitudes are useful, 



V2{0\diO)a^,j5Gap{vy)d{x)\7T^{q)) 



2(0|^^(0)^7^G„^(^;y)d(x)K"(g)) = Uq^ 



qx 



Paj(/)||(Qi)e" 



+ 



-V2{0\diO)^^j5Gap{vy)dix)\7r\q)) = Uq^ 



qx 



a,; e 



-iq(xa\+yva3) 



(12) 



where G^p = lea/SpaG^'^ , Vui = daida2da35{l - ai - a2 - as), g^p = ga/3 
The twist-3 and twist-4 light-cone distribution amplitudes can be parameterized as 



<t>37r{ai) 



360aia2a'3 



l + ^(7a3-3) 



30(5^ (ai — 02)03 



1 



+ 2e(l-2a3) 



<P\\i(^i) 



30<5^Q;i(l -03) 



a,; 



-120(5 010203 



+ 2e(l-2o3) 



- + e(l-3o3) 



(13) 



the nonperturbative parameters in the light-cone distribution amplitudes can be estimated 
with the QCD sum rules [SI |36]. 

After carrying out the operator product expansion near the light-cone, we obtain the 
following expression for the correlation function 11 H, 



n 



rribf; 



3lT 



27r2 

2 

'b 



dv I dt I Da. 



^0 



47r2 



dv I dt I Da. 











Mai)q -jp-k) [tvq • + 2(1 - t)(l - v)q ■ p] 
{m^ — [p — q{l — oi)]2} {m^ — [p — k — fOsg)^} 

q-{p-k) 

{m^ — [p — q{\ — oi)]2} {ml — {p — k — va^qY} 



[{l-2t) [2(l-?;),^^(o,)-(l-2^;)0||(o,)] +2(l-^;)0^(oi)-</'||(oi)} , 



(14) 



where = j^^- 

In order to suppress the contributions from the high resonances and continuum states, 
we perform n-th derivative with respect to the momentum {p — k)"^ to obtain stable n-th 
moment sum rules and Borel transform with respect to the momentum {p — q)'^ in the 
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i?-channel to obtain the Borel sum rules, then match with Eq.(lO), finally we obtain the 
sum rule for the nonfactorizable soft matrix element, 



2m 



Xcl 



[ ds I "^^^^^ / dt f da f 



dp 



^37r(l - a, a- (3,(3) 



tOc 

2(3 



1 1 
„ mf — mr, 

2 lo^-)_(x_t){\ 



ar^ mt - m: 



Xcl 



a 



a 



2 ptj pi pa 

M dt I da j 



d(3 



OLc. 



87r2 
(1 - 2t) 

arn^ + rn^^^ (1 — a) — 



2(l-^)0±-.^|| + 



2(l-^)0±-(l-^) 



exp 



aM2 



{1 — a,a — P, P) 



s + Q2 I a/? ' 



(15) 



where 



= - l-Wl 



4m2 



4m2 



On = 



t(l - t)s - ml 

t(l-t)(s-P2) 



— m 



(16) 



Xcl 



and Qo = ^m^^. In calculation, is chosen to be large space-like squared momentum 
(|P^| ~ ml) in order to stay far away from the hadronic thresholds, the value of cKc is a 
small positive quantity but not always small enough to be safely neglected, we perform 
the following approximation for the P integral. 



pa 
J ac 



dpG{s, X, a, P) 



J 



dpG(s, X, a, P) , 



(17) 



here G is an abbreviation for the integral functions and can be written as 



G{s,x,a,P) = A037r(l — a, a — /9, /3) + (1 — a, a — /3) + C(/)_l(1 — a, a — /?, /3) 
+L>(/)||(1 - a,a- P,P) +E(()±{1 - a,a ~ P,P) , 

A, B,C, D, E are formal notations. We can expand the light-cone distribution amplitudes 
(/•Stt) <;^|h (t>±-, (t^w and (t)±_ in terms of Taylor series of P, for example. 



b^T,{l- a,a- P,P) 



d 

(?!)37r(l -a,a- P, /?) 1/3=0 + ^'^37r(l - a, Q: - /?, P)\p=oP 
1 d'^ 

- - P)y=oP^ + ■ ■ ■ , (18) 



and continue into the timelike region analytically, = m^, then complete the integral 
J^" dpG{s,x,a,P). This procedure ensures the disappearance of the unphysical momen- 
tum k from the ground state contribution and enables the extraction of the physical matrix 
element due to the simultaneous conditions, P^ = m\ and {p — q)^ = m\. 
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The explicit expression for the physical hadronic matrix element is lengthy due to the 
re-summation of all the Taylor series of /3, here we show only the leading terms explicitly, 



An 



/ ds < — / dt da 



dp 



</)3^(l -a,a- P,P) 



ta 

Yp 



c / 2 

— s 



87r2 

(1 - 2t; 



1 

dt I da 



2(1 



«0 
P 



dp 



a 



ar 



(l-t)(l 



p 



a 



2(l--^)(/)±-0|| + 



(1 



2ac 



(1 — a, a — /?, /3) 



am. 



exp 



+ ml^^jl - a) - ml \ f m^^, 

aM2 M s + Ql 



n+l 



2m 



Xcl 



ds 



\ 47r2 



1 

aP 

1 

/ da 



dp 



^3^(1 - a,a - /3,/3) |/3=o 



toe . 2 



(1-0(1 



P 



^ 87r2 
(1 - 2t) 



if 



dt 



a 



da 



dp 







2(1 -^)0x 



1 + 



2(l-^)0.-(l-^ 



(1 - a,a - /?,/3) 



am. 



exp 



+ m2^^(l-a)-m2 \ / m2^^ +Q2 ' 
aM2 M S + Q2 



n+l 



a/3 



+ 



(19) 



In performing the P integral , we need only the values of the light-cone distribution 
amplitudes (/>3,r; 0|| , (t^\\ , 0± and their derivations at zero momentum fraction i.e. P = 0, 
there are no problems with negative partons (quarks and gluons) momentum fractions. 
The analytical continuation of to its positive value ends up with an unavoidable theo- 
retical uncertainty, if only a few terms of the Taylor series are taken. However, with the 
re-summation to all orders of P in Eq.(19), the assumption of quark-hadron duality is still 
applicable in the case of heavy meson final states. 



4 Numerical results 

The input parameters are taken as Vcd = —0.230, Vcb = 41.2 x 10~^, Vtb = 1.0, Vtd = 
8.1 X 10-^ m^^i = 3.511 GeV, tub = 5.2795 GeV [37J, fs = 0.18 GeV, sb = 35±2GeV2 
[32], /^^^ = 0.335GeV [33], sq = (16.0 ± 0.5) GeV^ [30j, U = O.lSGeV, mt = (4.7 ± 
0.1) GeV, nic = (1.35 ± 0.05) GeV, fs^ = (0.45 ± 0.15) x 10-2 GeV^, 003 = -1.5 ± 0.7, 
52 = (0.18±0.06) GeV^, e = f (0.2±0.1) [MUSSlEe] at the energy scale about /x = 1 GeV. 
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The parameters n and must be carefully chosen to warrant the high resonances and 
continuum states to be suppressed and obtain a reliable perturbative QCD calculation. 
The stable region for the Borel parameter is found in the interval = 11 it 2 GeV^, 
which is known from the B channel QCD sum rules [32]. In the charmonium channels, we 
usually perform n-th derivative and take n-th moment sum rules to satisfy the stability 
criteria 00], the ideal interval is n = 4 — 6. The parameter ^ is usually allowed to take 
values larger than 1 for the P-wave charmonia, we observe the best interval is ^ = 1.8 — 2.6. 

The nonfactorizable soft contributions come from the three-particle twist-3 and twist-4 
vr'^ light-cone distribution amplitudes, however, present knowledge about those distribution 
amplitudes is rather poor. The uncertainties of the nonperturbative parameters /stt, wa, 
5'^ and e are large, about (33 — 50)%, the nonfactorizable soft contributions An can be 
simplified into the following form, 

An = yif?,iv + y2hn^z + y?,5'^ + yA5'^e, (20) 

where the yi are numerical coefficients not shown explicitly. The uncertainties origin from 
the nonperturbative parameters /sjr, ws, 5^ and £ are rather large, even out of control. We 
take the same treatment as Refs. [201 [211 [221 [231 [26], and neglect the corresponding uncer- 
tainties, it weakens the predictive ability. The uncertainties origin from other parameters 
(sb, sq; "T-b; "^c, ^tc) are less than 22%, we take into account them with the formula 

A = fjr") (^i where the / denote the nonfactorizable soft contributions. 



the Xi denote the input parameters. 

Taking into account the next-to-leading order Wilson coefficients calculated in the naive 
dimensional regularization scheme p] for ^ = mft(mb) = 4.40 GeV and = 225 MeV, 

Ci(mF(mb)) = 1.082, C2(m^(mb)) = -0.185 , (73(m^(m,,)) = 0.014 , 
C4(m^(mf,)) = -0.035 , C5(m^(m;,)) = 0.009 , Ce{mE{mb)) = -0.041 , (21) 

finally we obtain the numerical ratio R between the contributions from the nonfactorizable 
and factorizable hadronic matrix elements, which is shown in Fig.l, 

p V2Ci(/i)(Xcl(p)^°(-g)|0(O)|g(p - q)) +0.20 

= r = U.b8_o,i6 • (22j 

The nonfactorizable soft contributions are rather large and they must be included in 
analyzing the branching fraction. The total amplitudes lead to the branching fraction, 

BriB^ ^ XciA = l-4lll]:i X 10-5 , (23) 
which is in agreement with the experimental data (1.12 it 0.25 it 0.12) x 10"^ marginally 

[a- 

The factorizable contributions of the operator ©2(0) can be analyzed in the same way 
with the correlation function Hpip, q, k), 



Hp = ^2 y d^^e-(^'-^)-y dV^^"')^(0|r{J,(y)O2(0)J5(x)}|7r0(g)) 

= I d^xe-'^P-'^^^0\T[d{0)rb{0)Mx)]\7r\q)) J d'ye'^P^'^^y {0\T[Jp{y)J,{0)]\0) 
oc {p-k),f{{p-kf)xfl^ + .... (24) 
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Figure 1: The ratio R between the contributions from the nonfactorizable and factorizable 
hadronic matrix elements. 

The contributions of the soft gluons can be absorbed into the i? — > vr" form-factor 
f {{p — k)'^) or the decay constant /xcd ^-nd differ from the sum rules for the operator 
0(0) greatly, where the contributions of the soft gluons are nonfactorizable, so we are free 
of double-counting. 

5 Conclusion 

In this article, we calculate the contributions from the nonfactorizable soft hadronic matrix 
element (xciip)'^^ {—q)\0 {0)\B (p — q)) to the decay — > Xci"^^ with the light-cone QCD 
sum rules. The numerical results show that its contribution is rather large and should 
not be neglected. The total amplitudes lead to a branching fraction which is in agreement 
with the experimental data marginally. 
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